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Some models describing unconventional fractional quantum Hall states predict quasiparticles that 
obey nonabelian quantum statistics. The most prominent example is the Moore-Read model for the 
v — 5/2 state, in which the ground state is a superconductor of composite fermions, the charged 
excitations are vortices in that superconductor, and the nonabelian statistics is closely linked to 
the degeneracy of the ground state in the presence of vortices. In this paper we develop a physical 
picture of the nonabelian statistics of these vortices. Considering first the positions of the vortices 
as fixed, we define a set of single-particle states at and near the core of each vortex, and employ 
general properties of the corresponding Bogolubov-deGennes equations to write the ground states 
in the Fock space defined by these single-particle states. We find all ground states to be entangled 
superpositions of all possible occupations of the single-particle states near the vortex cores, in which 
the probability for all occupations is equal, and the relative phases vary from one ground state 
to another. Then, we examine the evolution of the ground states as the positions of the vortices 
are braided. We find that as vortices move, they accumulate a geometric phase that depends on 
the occupations of the single-particle states near the cores of other vortices. Thus, braiding of 
vortices changes the relative phase between different components of a superposition, in which the 
occupations of these states differ, and hence transform the system from one ground state to another. 
These transformations, that emanate from the quantum entanglement of the occupations of single- 
particle states and from the dependence of the geometric phase on these occupations, are the source 
of the nonabelian statistics. Finally, by exploring a "self-similar" form of the many-body wave 
functions of the various ground states, we show the equivalence of our picture, in which vortex 
braiding leads to a change in the relative phase of components in a superposition, and pictures 
derived previously, in which vortex braiding seemingly affects the occupations of states in the cores 
of the vortices. 

PACS numbers: 73.43.-f, 74.90.+n, 71.10.Pm 



I. INTRODUCTION 

The experimental discoverji of the fractional quantum Hall effect (FQHE) led to intriguing theoretical observations 
regarding the elementary excitations (quasiparticles) of a two-dimensional electron system at a fractional Landau-level 
filling factor v. Very soon after the experimental discovery, Laughlin^ realized that the quasiparticles at filling factors 
v = l/(2p+ 1) (with p an integer) carry a fractional charge e* = ±e/ (2p+ 1) (for brevity, we use the term quasiparticles 
to refer also to quasiholes). Following that observation, Halperin^ showed that the hierarchy of observed FQHE states, 
at v = p/q (with q an odd integer), points to the fractional statistics of the quasiparticles and quasiholes, of the type 
that was previously studied by Wilczekii This observation was further clarified by Arovas, Schrieffer, and Wilczck. 5 
When a system contains two quasiparticles, and the positions of these quasiparticles are adiabatically interchanged, 
the state of the system acquires a geometric Berry phase. This phase, which is tt for fermions and 2tt for bosons, 
becomes a fraction of tt for FQHE quasiparticles. 

The experimental discovery^ of the even-denominator FQHE state v — 5/2 triggered the introduction of yet another 
novel concept with regard to the statistics of the elementary excitations. Employing conformal field theory to study 
the v = 5/2 FQHE, Moore and ReadZ discovered that if this state is well described by the Pfaffian wave function, 
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as numerical investigations seem to confirm^ the elementary excitations obey nonabelian statistics. The state of the 
system after a series of quasiparticle interchanges then depends on the order in which these interchanges are carried 
out. By using exact eigenstates of a model Hamiltonian^ Nayak and WilczekiS. subsequently showed that the ground 
state of the configuration in which 2N quasiholes are inserted at fixed positions is 2 JV -fold degenerate, and that the 
quasiparticles realize a 2 Ar ~ 1 -dimensional spinor braiding statistics. Following earlier observations that related the 
Pfaffian state to p— wave Cooper-pairing, Read and Greenii described this state as a p-wave BCS superconductor of 
composite fermions, and conjectured that the nonabelian statistics of its quasiparticles results from the zero-energy 
modes associated with vortices in this superconductor. 

This conjecture was further studied by Ivanov^ who mapped out the relation between the exchange of quasiparticles 
and the unitary transformation carried out on the Hilbert space of the ground states. Ivanov explicitly derived 
these unitary transformations, showed that they are indeed nonabelian, and confirmed that they are identical to the 
transformations derived earlier by Nayak and Wilczek 10 using conformal field theory. 

While these two derivations of the unitary transformations associated with vortex interchange may be easily gener- 
alized to calculate the transformations associated with other braidings of vortices, they do not provide a clear physical 
picture of the nonabelian statistics. This is exemplified in the following observation: using the method of Ref. Il2t it 
is easy to show that when the system is initially in a ground state |gs a ), and vortex j encircles vortex j + 1, then, 
under the assumption that no tunneling takes place between the vortex cores, the final state of the system is again a 
ground state, given by 

('•,< " + cje"* ') + ct^e"*"" 1 ) |gs Q ) , (1) 

where the operators , annihilate (create) a particle localized very close to the cores of the jth and (j + l)th 
vortex, respectively, and Clj is a phase defined in the next section. Eq. seemingly implies that the motion of the 
jth vortex around the (j + l)th vortex affects the occupations of states very close to the cores of the two vortices. 
This is in contrast, however, to the derivation leading to Eq. (I}, which explicitly assumes that vortices are kept far 
enough from one another so that tunneling between vortex cores may be disregarded. 

In this work we study the ground states |gs Q ) in an attempt to give a physical picture of the effect of braidings in 
the positions of vortices. We show that two ingredients are essential for the nonabelian statistics of the vortices. The 
first is the quantum entanglement of the occupation of states near the cores of distant vortices. The second ingredient 
is familiar from (abelian) fractional statistics: the geometric phase accumulated by a vortex traversing a closed loop. 

Within the Chern-Simons composite-boson theory^ the abelian fractional statistics of the v = 1/m states is ex- 
plained by mapping the ground state of the electronic system to a superfluid of composite bosons, and the quasiparticle 
excitations to vortices in that superfluid. Due to the coupling of the vortex to a Chern-Simons gauge field, the de- 
pletion of bosons at the vortex core is quantized to a fraction 1/m of a fluid particle. Thus the charge carried by the 
vortex is also fractional. The quantum statistics is related to the geometric phase accumulated by a vortex traversing 
a close trajectory. Roughly speaking, the vortex accumulates a phase of 27r per fluid particle which it encircles. When 
another vortex with its fractional charge is introduced to the encircled area, this phase changes by a fraction of 2-k. 
Upon adapting the argument to interchanging of vortices, one finds that this fraction of 2-k translates into fractional 
statistics. 

Similarly, the Moore- Read theory of the v = 5/2 state describes it also as a superfluid, with the quasiparticles being 
vortices in that superfluid. However, the "effective bosons" forming the superfluid are Cooper pairs of composite 
fermions. Consequently, the superfluid has excitation modes associated with the breaking of Cooper-pairs. In the 
presence of vortices, a Cooper-pair may be broken such that one or two of its constituents are localized in the cores 
of vortices. For p— wave superconductors, the existence of zero-energy intra-vortex modes leads, first, to a multitude 
of ground states, and, second, to a particle-hole symmetric occupation of the vortex cores in all ground states. When 
represented in occupation-number basis, a ground state is a superposition which has equal probability for the vortex 
core being empty or occupied by one fermion. 

When a vortex traverses a trajectory that encircles another vortex, the phase it accumulates depends again on the 
number of fluid particles which it encircles. Since a fluid particle is, in this case, a Cooper pair, the occupation of a 
vortex core by a fermion, half a pair, leads to an accumulation of a phase of ir relative to the case when the core is 
empty. And since the ground state is a superposition with equal weights for the two possibilities, the relative phase 
of 7r introduced by the encircling might in this case transform the system from one ground state to another. 

This qualitative picture is made more precise in this paper. Our analysis revolves around the definition of a set of 
single-particle states localized at or near the 2N vortex cores. We start by defining the "core states," a set of 27V states 
each of which is localized at a specific vortex core. We find that in all possible ground states, the occupation of these 
single-particle states near one vortex is entangled with the occupation of single-particle states near all other vortices. 
We prove that any many-body state in which these occupations are disentangled is necessarily an excited state. We 
show that the evolution of the ground state as positions of vortices are braided indeed follows the picture outlined 
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above, and discuss both the case where vortices encircle one another and the case where they interchange positions. In 
making this picture of nonabelian statistics more precise, wc define, starting from each core state, further orthogonal 
single-particle states ("near-core states") which are localized near a vortex core. The occupations of these additional 
single-particle states in the many-body ground states are also particle-hole symmetric and entangled between different 
vortices. In fact, we reveal a "self-similar" structure of the many-body wave function with respect to the occupation of 
these single-particle states which leads us to express the relation between our picture of nonabelian statistics and the 
known representations of vortex braiding in the space of ground statesifi*i£ in terms of compact operator identities. 

The structure of this paper is as follows: We begin in Sec. [H] with a review of the description of the Pfaffian 
state as a p— wave superconductor, with vortices as quasiparticles. In Sec. IIIII we start with the definition of the 
single-particle states by introducing the "core states." In Sec. II VI we explore the roles of quantum entanglement and 
geometric phases in the evolution of these superpositions when vortex positions are braided. The "near-core states" 
are introduced in Sec. E] The occupations of these "near-core states" in the many-body ground states is worked out 
in Sec. I VII revealing the "self-similar" structure of the wave functions. We conclude in Sec. IVIII Some details are 
relegated to appendices. 



II. SOLUTIONS OF BOGOLUBOV-DEGENNES EQUATIONS - REVIEW 

The Pfaffian trial wave function for the quantized Hall state at Landau level filling factor v = 5/2 was first introduced 
by Moore and Read- as the first-quantized wave function 



W*i>**"0 = Pf — r I[(*-*i) 2 II e > ( 2 ) 

\ Z i Z 3ji<3 3 



where t is the magnetic length and z% — x% + %y% is the complex coordinate of the zth particle. For p particles, the 
Pfaffian in Eq. takes the explicit form 



Pf = n W, -4 <^ \ , (3) 

\z t - ZjJ 2Pi z {p/2)\ { zi - z 2 z 3 - z A Zp-i - z p J 

where A is the antisymmetrization operator. It is instructive to view the Pfaffian appearing in the wave function in 
Eq. J2J) as the real-space BCS wave function of composite fermions for a fixed number of particles>ii*i^ According 
to the associated pair wave function g{z) — 1/z, the pairing is of spinless (or spin-polarized) composite fermions 
in the I = — 1 angular-momentum (p-wave) channel. The Pfaffian corresponds to a weakly-paired superconductor, 
which for a two-dimensional p-wave superconductor is topologically distinct from the strongly-paired phased The 
charged excitations of the quantum-Hall system are, in this description, the half- flux-quantum (h/2e) vortices of the 
superconductor. 

As is often the case, the use of a first-quantization formulation is mathematically involved, and makes the physical 
picture difficult to read. The identification of the Pfaffian as a complex p-wave BCS state of composite fermions 
subsequently led Read and Greenii to introduce a second-quantization formulation which paved the way for a clearer 
physical picture. Their starting point is the BCS mean field Hamiltonian 

H = J drV> t ( r )M'(r) + 5 / dr dr' {D*(r, r')ijj(r')ijj(r) + D(r, r')^ (r)^ (r')} (4) 

with the single-particle term ho and the complex p-wave pairing function 

D(r, r') = A (l±*) (id x , - d y ,)5(r - r'). (5) 



Read and Greenii retain only the potential part of ho by setting ho = — /ti(r) and argue that this is sufficient for 
studying the topological properties of the Pfaffian state, such as the statistics of its quasiparticles. In the presence of 2 N 
vortices pinned at positions Hi, the gap function takes the form A(r) = | A(r)| exp[i%(r)] with x( r ) = Si=i arg(r— Rj). 
In the vicinity of vortex k, the phase x( r ) can be approximated by x( v ) — ar g( r ~ Rfc) + ^fc with flk — J2i^k ar g(Rfc — 
Ri). 

The fermionic excitations of superconductors are described by the Bogolubov-deGennes (BdG) equations 
F (u(r)\_( -fi(r) i{A(r),d x +id y }\ fu(r)\ 

*\v{t)J ~ \i{A*(r),d x -id y } M (r) )\ «(r) J ' W 
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For two-dimensional complex p-wave superconductors, solutions of non-zero energy should be distinguished from 
those of zero energy. We denote the non-zero energy solutions by (u E (r), v E (r)), with u E (r) = v*_ E (r). In 
second quantization, positive-energy solutions are associated with annihilation operators of BCS quasiparticles 
Tfi = / dr [ue(t)iP(t) + ve(t)i/)^ (r)\ , while negative-energy solutions are associated with creation operators of the 

same quasiparticles, T- E = F E . The zero-energy solutions (ui(r),Vi(r)) are localized at the vortex cores. For well- 
separated vortices, there is one such solution per vortex. With a choice of phase the advantages of which become 
clear below, these zero-energy solutions take the form 

Ui (r)=v*(r) = ^=w^ (r)e^. (7) 

Here, the index i = 1, . . . , 2N numbers the vortices and the functions u>^(r) are normalized wave functions localized 

near the core of the ith vortex. When the vortices are well separated, the functions (r) are mutually orthogonal. 
In Sec. [V] below, we iteratively define additional single-particle states localized in the vicinity of the vortex cores. 
These states will be denoted by (r) , with the subscript labelling the vortex and the superscript enumerating the 
states near each vortex. 

The zero-energy eigenstates of the BdG equations correspond to the Bogolubov operators 
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(8) 



where we introduced the operators Cj = J drwj°' ) (r)-0(r) which create particles in the vortex-core states Wj(r). 

Evidently, 7J = jj so that the operators associated with the zero-energy solutions are Majorana fermions. 

The existence of the zero-energy solutions leads to a degeneracy of the ground state. Enumeration of the ground 
states is customarily done by combining the Majorana operators 7$ in pairs and defining (ordinary) fermionic creation 
and annihilation operators 

a lj = \ (723-1 + *72i) (9) 

with j = 1, . . . , N. The ground states can now be written in the occupation- number basis corresponding to these 
fermionic operators. A ground state |m) is then labelled by the occupation numbers m = (7712,7714, . . . ,TO2At) with 

4 i a 2j -|m) =m 2 j|m) (10) 

and 

|m) = (a\ N r™ . . . (4r> 2 r 2 |m = 0) (11) 

leading to a 2^-fold degeneracy of the ground state. 

The BCS Hamiltonian Eq. Q is diagonal when written in terms of the quasiparticle operators T E and 7$ (Bogolubov 
transformation). The ground state is determined by the conditions |gs) = for all E > 0. For a uniform 
superconductor in the absence of vortices (i.e., for space-independent A and /i), these equations lead to the celebrated 

BCS wave function |BCS) = Yl^ ( u k + v ^ c L c ^--kj l vac )i where the prime indicates that pairs of momenta (k, — k) 
should be counted only once, and |vac) is the state with no particles. In this case, the choice of a plane-wave basis 
for the single-particle states is natural. 

In the presence of vortices, A and /1 are space dependent, and there is no obvious single-particle basis for the 
description of the ground states. A proper choice of such a basis turns out to be helpful in our discussion of nonabelian 
statistics. 



III. CORE STATES 

A natural starting point for a single-particle basis are the 2N states w\°\r). These single-particle states, which we 
approximate to be orthogonal due to the large distance between the vortices, are associated with a 2 2Ar -dimensional 
subspace of the many-particle states. The remaining (infinitely-many) single-particle basis states remain unspecified 
throughout this section and are partially defined in Sec. We will refer to the first 2N single-particle states as the 
"vortex-core states," and to the remaining single-particle states as the "other states." 
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Many-particle states in Fock space can now be expanded in terms of occupation numbers of these single-particle 
states. A corresponding basis state of the Fock space is then written as 

vortex core other (12) 

where (1) denotes an empty (occupied) state. The first factor has 2N digits and designates the occupations of the 
vortex-core states. We enumerate these states by |t) with t — 1, . . . , 2 2N . The second factor in Eq. I|12(l designates 
the occupations of the "other states". Below we find that for every possible ground state |gs a ) the probability to find 
the core states in an occupation |t) is equal to 1/2 2N , and is independent of t. A state in which this occupation 
depends on t is necessarily an excited state. 

Although the operators a 2 j and aL act on the occupations of the core states only, the quantum numbers m do not 
fully label the 2 2N — dimensional space of states |t). Specifically, there is no direct relation between the occupation 
numbers m and the occupation of the single-particle states Wj(r). In order to explore the structure of the ground 
states in terms of the states |r), we now introduce another set of quantum numbers associated with a second set of 
N Majorana operators, 



•v, ^(<y " c}< (13) 

whose associated BdG spinors (iw^°\r)ei nj , —i[w^°\r)]*e~^ sij ) are orthogonal by construction to the zero-energy 
solutions of the BdG equations. These operators obviously also act only on the occupations of the vortex-core states. 
However, when expanded over the complete set of BdG quasiparticle operators, they involve only non-zero energy 
quasiparticles so that 



£[c°r B + (w) 



_E>0 



with coefficients C 3 E — i\[2 J dr u^(r)Wj (r)e l0j / 2 . Upon pairing the vortices, these Majorana operators can again 
be combined to obtain (ordinary) fermionic operators^ 



^ 2j = ^=(iX 2j _ 1+ X 2j ). (15) 

We label the occupation numbers of these fermions by x 2 , X4,, . . . , x 2 n, and introduce x = (x 2l X4, . . . , x 2 n). 
We now form a basis of the 2 2Ar -dimensional Fock subspace of the vortex-core states by defining 

|m, x) = {a\ N r^ . . . (a\r*(4r* ifll)** GSt)** . . . tf\ N )**» |m = 0, x = 0) (16) 



where |m = 0, x = 0) is the state annihilated by all the a 2 j's and f3 2 /s. We obviously have (m, x|m', x') = <5 mm '<5 xx < 
In terms of the occupations of the vortex-core states, the states (|16f) take the explicit formic 



N 

|m,x) = wIl{[ 1 + ! (- 1 ) l2,e ^ (!!!i " +!!!,> 4 J -i4 J ]^ 

+ [e-^-c^ + i (-IJ^e-^ctJ <i„,. .,..,} |vac) . (17) 



■m 2 j,x 2 j 



The sign factor s m , x = Yli=2 l) X2,( - m2r+X2r ^ arises due to the different operator orderings in Eq. (fTCjl and in 

the product in Eq. (|T7|) . 

A ground state labelled by m is then a superposition of the form 

|m) = £|m,x)|A x >. (18) 

X 

It is important to note that the states \A X ) are independent of the particular ground state |m). Arbitrary ground 
states |gs a ) can be written as linear superpositions of the states |m). Here and below we use the notation \A) to 
denote states in the Fock subspace corresponding to the unspecified "other states" in the single-particle basis. There 
are 2 N components in the superposition (|18f) . one for every value of x. The states |^4 X ) should be determined by the 
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requirement that the ground states are annihilated by all positive-energy annihilation operators r^. Although we do 
not know the complete set of operators , we can now show that 

(A X \A*) = ^5^. (19) 

To see that, we first note that since the operators Xj are composed of finite-energy quasiparticle operators only, 
the matrix element for any odd number of X operators between any two ground states must vanish. Thus, 

<gs Q \X n | g S/3 ) = <gs Q \X h X j2 X h | g S/3 ) = . . . = (20) 

for arbitrary indices a,fi,jx,... Second, the matrix elements of a product of two different operators Xj between states 
in the ground-state manifold are 



< g s Q \x h x h \ gS0 ) = 8 aP ]T c£cr 

E>0 

= 2S aP [ drdr' U ;j° ) (r)[u,g ) (r / )]*eK fi «-^) £ u%(v)u E (v'). (21) 

J E>0 

For sufficiently high energies the functions ue(y) are approximately plane waves so that J^e u *e( v ) u e{ t ') is a short- 
ranged function, presumably decaying exponentially with |r — r'|, even for non-uniform superconductors. Then, for 
well-separated vortices ji and jit the matrix elements in Eq. (|21|l approximately vanish. Put in different words, the 
operation of the operators Xi is spatially localized around vortex i, and thus two such operators operating near two 
distant vortices generate orthogonal excitations. Similarly, the matrix element of any other even number of different 
X operators, taken with respect to any two ground states, vanishes as well. In Sec. [V] we will also give a direct 
algebraic proof of this result which relies on a relation of the Xi to operators annihilating the ground states. 

The conditions in Eqs. (|20|l and (|21fl imply some general conclusions regarding the states \A X ), which we first 
present for the case of two vortices. There are four states |r), labelled by 1 00) , 1 01) , 1 10) , The phases encoding 

the vortex positions are related by 0,2 = ^1 + tt. The two ground states then take the explicit form 

|m = 0) = (|00)-e- 40l |H»l^o)+e-^ l/2 (|10)-|01))|A 1 ) (22) 
|m=l) = (|00)+e- 4Ol |ll»l^i>+e^ l/2 (|10) + |01))|A ). (23) 

The condition (gs a \Xi\gs a ) = implies that (^ol^i) = 0. In addition, (gs a |XiA"2|gs Q ) = imposes (Ao|Ao) = 
= 1/2. Thus, while we cannot find the complete wave functions of the ground states without a full solution 
of the BdG equations, our procedure leads to the conclusion that the two ground states are incoherent superpositions 
of the states (|00) ± e~ lQ-1 |11)) with the states e~ %rtl ^ 2 ( 1 10) ± 1 01) ) , with equal weights to both components. There 
is an equal probability 1/4 for all four possible charge arrangements 1 00) , 1 01) , 1 10) , |11). The parity of the particle 
number differs between the two basis vectors of the ground-state manifold. However, this difference in parity does 
not originate from the occupations of the two vortex-core states. A local measurement of one of the two vortex-core 
states cannot distinguish between the two ground states. 

Generalizing to 2N vortices, it is easy to see that the requirement that the expectation value of products of 
X operators vanishes leads to Eq. Q19[l. Thus, all basis functions of the ground-state manifold are an incoherent 
superposition of 2 N terms, of equal weight. Each of these terms is, by itself, a coherent superposition of 2 2JV ~ 1 
possible occupations of the core states, which constitute all possible occupations of a given parity. This observation 
clarifies why there are 2 N ground states, rather than 2 2N . Within the Fock space of the core states there are 2 N 
eigenvectors for each eigenvalue m, but a ground state is generated by their incoherent superposition, mixing all of 
them with equal weights. 

The wave functions in Eq. (|18fl give an explicit description of the occupations of the core states. The operator 
7V(uij- ') = (1/2) ('ji + iXi) ('ji — iXi) counts the number of particles in the ith core state. It is easy to see that 



(m'|iVK (0) ) |m) = (m'| [n(w^)\ ' |m) = l -8 mm ,. (24) 



Thus, in all possible ground states (including arbitrary superpositions of the states |m)) the occupation of each core 
state is particle-hole symmetric, with a probability of one half for being empty or occupied. 
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IV. GEOMETRIC PHASES AND QUANTUM ENTANGLEMENT IN THE EVOLUTION OF CORE 

STATES UNDER VORTEX BRAIDING 

When vortices are adiabatically moved along closed trajectories, ending with braiding of their positions, the ground 
state may evolve in time away from the initial state, and the final ground state may thus be different from the initial 
one. There are, in principle, two contributions to this transformation of the ground stated The first originates from 
its explicit multiply-valued dependence on the phases £li (explicit monodromy). The second is through a nonabelian 
geometric vector potential that gives rise to a nonabelian Berry phase. As is always the case with geometric phases, 
only the sum of the explicit monodromy and the Berry phase is observable and one can split this sum between both 
contributions in any way desired by choice of appropriate phase factors. We now show that the phase choice we made 
in Eq. Q makes the second contribution vanish, and proceed to calculate the first contribution. 

We need to prove that the geometric vector potential-Si 

Im(m|V R Jm') (25) 

vanishes for all m,m',Ri. The states |m) depend on Rj through the phases fli and through the functions u>|°^(r). 
Since the w[ (r) are real (up to some trivial global phase factor) , their derivative does not contribute in Eq. i|25|l , 
and we may write Vr ; = 5Zj(Vn. i f2 i j) and compute the matrix element 



<m|-^-|m') = «W$> X | 75 ^-|A X )+ * £< m ,x|-^|m',x). (26) 



The first term on the right-hand side is diagonal in m, m' and otherwise independent of m. It therefore leads only to 
abelian phase factors. Using the explicit states in Eq. (|17|l one finds that also the diagonal elements m = m' of the 
second term are independent of m. The only contribution to the nonabelian part of the Berry phase can therefore 
arise due to the off-diagonal contributions m ^ m' of the second term. However, these vanish as one may verify by 
using the explicit states in Eq. 1)170. Thus, when the wave functions are written as in Eqs. H17f) and H±8|l the only phase 
factors that lead to non-trivial unitary transformations arise from the explicit dependence of these wave functions on 
the phases Qj. 

More explicitly, Eq. I|17|l shows that the part of the wave function in which the core state i is occupied has an 
amplitude proportional to e - if W 2 ; while the part in which this state is empty does not depend on flj. Since the ground 
states involve superpositions of empty and occupied core states, when the phase Oj accumulates a 2tt shift, a relative 
minus sign is introduced between the components of the superpositions in which core state i is empty and occupied, 
and the ground state does not necessarily come back to itself. The source of the evolution from one ground state to 
another is then in the phases between different components of the wave function, and not in any change of occupation 
of states. 

It is instructive to examine in detail the case of the ith vortex encircling the jth vortex, for which 

Qi -> + 27r 

flj -> flj + 2tt. (27) 

This change of phase affects both the states |m,x) and the states \A X ), according to 

|m,x) -> (2i-f i X i )(2i 7j X j )\m, X ) (28) 
\A X ) -> 4 l^x> (29) 

where By is an operator that acts on the "other" , non-core, states, only. The final state 
2 x (2i7jXj)(2i7jJCj) |m, x)By \A X ) must be a ground state. However, the operators X i: Xj generate excitations 
above the ground states. The operator Bij must then be an operator that annihilates these excitations. More pre- 
cisely, since the states \A X ) do not depend on m, the operator XiXjBij must be a c-number within the subspace of 
ground states. In fact, for the norm of the ground state to be conserved during the braiding of vortices, the magnitude 
of this c-number must be unity, i.e., 

2X l X j B lj = e^" (30) 



with (f> a being an abelian phase. In this way, we recover the known unitary transformation for winding of two vortices, 
^■lilj given in Eq. [see also App.Ia"]. This shows that despite the appearance of Eq. QJ, this transformation does 
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not involve any changes of core-state occupations, but only changes in relative phases between the components of a 
superposition, each of which has different core-state occupations. 

In the subsequent sections, we will study the states \A X ) much more explicitly by introducing the "near-core" 
single-particle states. We will find that these states are structurally very similar to the ground states |m) and this 
allows us to give explicit expressions for the operators Bij (up to an abelian phase factor). These expressions do 
indeed satisfy Eq. l(3Tjl) . 

The above observations allow us to conclude that the ground states spanned by the basis in Eq. (|18|l are states 
in which the occupation of the core states at different vortices are fully entangled, in the following sense: There is no 
ground state of the system in which the parity of the number of particles is well defined, and in which the occupations 
of two subsets of core states are disentangled from one another. If such a state were to exist, its wave function could 
be written as 

C\C 2 C 3 |vac) , (31) 

where C\ is an operator that acts only on the core states belonging to the first subset, C 2 is an operator that acts only 
on core states belonging to the second subset, and C 3 is an operator that acts solely on the remaining states (core 
states or other states) not included in any of the two subsets. However, for the state in Eq. (|3 II) to have well-defined 
particle-number parity, the states created by C\ or C 2 must each have definite parity. By considering the effect of an 
encircling in which one vortex (say, the ?th vortex) from the first subset winds around another vortex of the second 
subset (say, the jth vortex), we find that this cannot happen. The unitary transformation 27^ corresponding to this 
transformation changes the parity of the particle number for both the first and the second subset, while we showed 
that the transformation of the ground state is a consequence of changes in phases only, rather than changes in core 
states occupation. Consequently, a state of the form in Eq. (|31[) cannot be a ground state. 

In some sense, the entanglement of the occupations of core states is signaled by Eq. Hl()(l . As seen from that 
equation, the ground states 72.7-1 |m) and j 2 j |m) differ from one another only by a phase factor, despite the fact 
that they are obtained from the ground state |m) by the application of two Majorana operators localized very far 
from one another. 

Interestingly, despite the entanglement, the occupations of different core states i and j are uncorrelated, 

(gs|iVK (0) )iv(^ 0) ) |gs) = (gs|7vK (0) ) |gs) (gs\N(wf) |gs) = ~. (32) 

This lack of correlations persists also to higher-order correlators. 

Furthermore, we can conclude that the total number of particles in the core states, counted by the operator 
Arorc = Ym=i ^( w i°^)t does not have a well-defined parity. Rather, it has equal probabilities for being even and odd, 
irrespective of what is the parity of the total number of particles in the ground state. To see that, we consider the 
particle-number parity operator exp(i7riV C orc)- Due to the lack of correlations between different vortices, 

2 N 

(gs| exp (i 7 rJV core )|gs) = [] <gs| exp [mN{wP )] |gs) = 0. (33) 
i=i 

Thus, the parity of the ground states cannot be determined by a measurement of the occupation of the core states 
alone. 

Despite the inherent quantum entanglement of the ground states, it appears impossible to formulate corresponding 
Bell inequalities. Each core state defines a two-dimensional Hilbert space, similar to a spin-^. However, the operators 
associated with these spaces at different cores, the fermionic operators jj and Xj, do not commute. Thus, their 
classical analogs are ill-defined^ 

What happens when two vortices are interchanged can be analyzed along similar lines to the case of one vortex 
encircling another. We now analyze the case of interchanging vortices from the same pair, 2i — 1 and 2i. Note that 
choosing the vortices from the same pair does not imply any loss of generality, both because for any interchange we 
can choose a pairing such that the two vortices are from the same pair, and because our considerations will eventually 
lead to an operator expression which is independent of the particular choice of pairing. 

As the adiabatic motion of vortices does not involve any tunneling of particles between vortex cores, an inter- 
change of the position of vortices interchanges the occupation of their core states. When both core states are 
occupied, the interchange is accompanied by a factor of (—1), since two fermions interchange positions. In addi- 
tion, one of the vortices necessarily crosses the cut line of the phase of the other, changing this phase by 2tt (cf. 
App. f^). The phase of the other vortex remains intact. Implementing these transformations in Eq. I|17|) . the term 
[1 +i (— l) X2i e~ l( - n2j - 1+n2j ^ 2 C2j_ic] !: j} remains unaffected by the interchange, while the relative sign between the two 
components changes in the term [e~ in2j - 1 ^ 2 cl ■_ 1 + i (— l) X23 'e~* n23 '/ 2 Cr, •]. 
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These transformations of the core-state wave functions are implemented by the operator e 71 ' 72 ^ 172 ^ I" 1 e nX2 i- lX2 i / 2 , 
acting on Eq. I|17|). By contrast, as found in Ref. ifl^ (reviewed in App. |A"I). the unitary transformation that enacts 
the vortex interchange on the states in Eq. Ijl8|l is e 71 " 72 ^ 172 ^ 2 . As in the case of encircling, the difference between 
these two transformations is an operator B that acts on the \A X ) part of the wave function in Eq. (TSJ. 

The states \A X ) together with the operators B that act on them when vortices wind and interchange are the subject 
of the following sections. We define a set of single-particle states adjacent to the core state for each vortex ("near-core 
states" ) , and show that the operator B affects the occupation of these states in exactly the same way as the operators 
27i7„- and e 7r72 ^ l72j / 2 affect the occupation of the core states for vortex winding and encircling, respectively. In fact, 
we find the corresponding operators B to have the same functional forms, but with the Majorana operators j{ and 
7j replaced by analogous Majorana operators associated with the near-core states. We explain this in terms of a 
"self-similar" structure of the ground-state wave functions. 



V. STATES NEAR THE CORE 

Our information on the states \A X ) introduced in Eq. (|18J) has so far relied on the conditions in Eqs. H20J) and Q21|l 
for matrix elements of products of the operators X^. These operators were constructed to have two main virtues: 
(i) They create and annihilate particles in the core states only, (ii) They create and annihilate only excitations 
(Bogolubov quasiparticles) with E > 0. We now take these operators as a starting point in a scheme by which we 
define the set of mutually orthogonal single-particle states w^(r), where the subscript i refers to the vortex number 
and the superscript k enumerates the iterations in our scheme. The new single- particle states remain localized near 
the vortex cores and provide insight into the nature of the states \A X ). 

(k) 

We get to these states by using the Aj's to construct a set of operators that, unlike the Majorana Aj's, 
annihilate the ground states, 

Y t (k) |m) = (34) 

for all i,k, m. Furthermore, these operators anti-commute {if , if } = and each operator if' creates and 

annihilates particles only in the two basis states w\ k (r) and (r). By virtue of the conditions in Eq. (|34|l . these 

operators specify the ground-state wave functions in the subspace spanned by the states loj (r) . 

Our scheme starts from the expansion of the operators Xi in Eq. (|14|l . We define the two sets of operators 

yf } = V^^C^Te (35) 

E>0 

Z P = * E ( C e T z - °e T e) ■ (36) 

E>0 

The operators if annihilate all the ground states since they are constructed from (positive-energy) annihilation 
operators only. The operators z\ are, by construction, Majorana operators, and can thus be written as 



z >" = 7i hh 



1 



w 4 (1) (r)^(r) +h.c. . (37) 



This defines a set of 2N single-particle states to^ (r). We prove in appendix iBl that for well separated vortices the 
w\ (r) are mutually orthogonal as well as orthogonal to the core states wf 3 (r). Thus, we can extend our single- 
particle basis by adding to it the 2N states u>-^ (r). Note that for conciseness of notation, we absorb the phase factors 

exp(if2j/2) into the definition of the states it?j (r) throughout this section. We comment below on how they should 
be reinstated. _ 

Since X t = (l/\/2)K (1) + Y^}, the condition 

if > ISO = (38) 

implies immediately that any combination of different operators Xi has zero matrix elements between ground states 
as required in Sec. EH Using the relation Y t (1) = (l/V2){X l - iZ^), we can write 

Y! = ~ / dr {( Wl (0) (r) - Wl (1) (r)) ^(r) - ([w^ (r)]* + K (1) (r)]*) ^ (r )} . (39 ) 
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Thus, these 2N operators affect the occupations of the states w^(r) and w^(r) only. 

To iterate this process it is helpful to summarize the steps leading to the definition of w^\r) through Y t ' , us- 
ing the concise spinor representation for the operators. We start with the operators 7, = A^{w^ (r), [w4°^(r)]*) 



in Eq. (JSJ) which act on the occupation of the wl 0> (r) only. In Eq. (|13|) we define the operators Xi 



-^(j- w \ ( r ), ~i[ w i ( r )]*) acting on the occupations of the same single-particle states (a z is a Pauli matrix). By con- 
struction, the spinors corresponding to Xi are orthogonal to those corresponding to 7,; , implying that the corresponding 
operators anti-commute. Then, in Eq. I|35|) . we extract from the operators Xi the parts that annihilate the ground 
state. The operators Yi, in turn, are written in Eqs. i|36|) - (|39[) as sums of hermitian operators X{ and anti-hermitian 



operators iZ^, and finally, we defined ^^(r) through the Majorana operators Z\ L) . In spinor representation, these 
last three steps can be recast as Z^ = —Sa z "fi with 



(i) 



S = s gnE 



U E 
VE 



UE 
VE 



(40) 



The operator S is a difference of two projection operators. The positive (negative) energy part of the sum projects 
to the subspace of positive (negative) energy BdG solutions. 



We now iterate this process to define states w\ k '(r) and a set of operators Y^ K> for which Y^ K> |gs Q ) = 0. This is 
achieved by generating the set of Majorana operators 



(A') 



Z 



(fc) 




[Sa z ] ' 7i 



(41) 



When writing out S explicitly in real space representation, according to Eq. Q40[l. it contains energy sums of the 
type that has been discussed following Eq. J2TJ. By the same arguments employed there, we can conclude that as 
a function of the two coordinates r and r', the operator S is short ranged, i.e., decays fast for large |r — r'|. Thus, 
(r) is also localized around the ith vortex although its distance from the vortex core increases with k. Since the 



construction assumes that states (r) localized around different vortices do not overlap, there exists an upper limit 
to the number of iterations. We denote the last iteration number by L. 



The functions w\ (r) are studied in more detail in App. [51 We find that for iterations k < L, the various w\ (r 
are mutually orthogonal, 



K K ) 



SiiS, 



ij Okk' 



(42) 



As long as the spatial extent of the states (r) is small compared to the distance between vortices, the dependence 

of the near-core states (r) on the phases f2j is analogous to that of the core states. Specifically, these phases can 
be made explicit by the replacement 



Wl (fe) (r) 



w 



(fe) 



(r)e 



iQi/2 



(43) 



This can be seen from Eqs. H40|) and (|41|) in combination with the observation that for r in the vicinity of vortex i, 
all finite-energy BdG solutions ue(t) and ve(i") depend on the vortex positions through the phase factor e in ^ 2 . 
We also define 



1 w. 



(fe) 
(fe) 



(r) 



(44) 



and the operators 



l.(l + S)a z Zt 1) 



(45) 
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which annihilate the ground states [cf. Eq. J33J]. Here, the operators annihilate particles in the states w[ k \r). 

This iterative construction enlarges our single-particle basis. Starting with the 2N core states w!- \r), we defined an 
additional L — 1 "near-core" states for every vortex. 

The occupations of the newly-defined L — 1 "near-core" states are all particle-hole symmetric. This can be seen by 
noting that the operator N(w^) = (1/2) (^Z^ + iX- k ^ (^Z^ — iX- k ^ counts the number of particles in the state 

w-^(r). Since 



1 + 

1 

- 1 



+ (^ (fc+1)t +^ (fc+1) ) ( 



y(*)t _ y W 



and since the operators Y^ k ^ and Y^ k+1 ^ correspond to orthogonal excitations, one finds 



m'| N{w {k) ) |m) = (m'| N(w 



|m) = - 8 mm , 



(46) 



(47) 



In the next section we solve Eq. (|34fl to extract the ground-state wave functions in the Fock subspace spanned by 
the near-core states '(r). 



VI. GROUND-STATE WAVE FUNCTIONS AND NONABELIAN STATISTICS 

Solving Eq. (|34|) . we now show that the structure we found in Sec. 11111 for the occupations of the core states repeats 
itself for the near-core states. This allows us to complete the argument that the effect of vortex braiding can be 
understood explicitly in terms of the geometric phases . 

The detailed form of the ground-state wave function depends on the convention for the order in which creation 
operators of various states act on the vacuum. Eq. H18fl positions all creation operators for the core states to the left 
of creation operators for other states, and Eq. ljTT)|) specifies the order in which core states are being filled. In that 

spirit, the creation operators for the first iteration of near-core states w^\r) are going to be positioned to the right 
of those of the core states, and so on with increasing number of the iteration. Within each iteration, we follow the 
convention defined in Eq. Ijl6(l for the core states. 

We first determine the occupations of the near-core states of the first iteration, (r) in the states |A X ). To this 
end, we define creation and annihilation operators from pairs of the operators zf^ and X\ ; , 



5 2j = -^(Z^-iZ^) (48) 
V2j = -^(-iX$_ 1+ X$). i 4') ) 



V2 

J_ f _ ; yU . r U 

Then we can express the states \A X ) in terms of the states 

im,x)« = (4 N r™ ■ ■ ■ (6ir*(4y n *(4r*(r,ir* ■ ■ ■ (4 N ) X2N \™ = o,x = o>« (50) 

where the superscript on the state indexes the iteration. The conditions (Y^j-i ^ r )l m ) = f° r ^ nc g roun d states 
can be rewritten as [cf. Eq. (|43|) ] 

G^y- 4,01m) = (51) 

(^.+<5 2j )|m) = 0. (52) 
Applying these conditions to the ground states, we find 

|A x )=]>>,x')«|A x ,}« (53) 

x' 

in terms of further states l-Ax*)^ which contain the occupations of the remaining near-core states as well as the other 
states. As expected, Eq. (|53|) satisfies Eq. (|19ll . Iterating this procedure, we get 

|A x )«=£|x,x')( 2 Vx'> (2) , (54) 
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and so on, as long as near-core states from different vortices do not overlap. Consequently, we can write the ground- 
state wave functions as 

|m) = £ |m,x(°))(°) £ |x(°),x (1) > (1) E l* (1 U (2) > (2) ■ ■ ■ £ Ix^x^)^)^,)^). (55) 

x (0) x (l) x (2) X (L) 

This structure of the ground states is helpful for analyzing the effect of vortex braiding. 

As discussed in Sec. IIVI the effect of encircling vortex i by vortex j on the ground states |m) is given by the 
unitary transformation 27^ [Eq. QJ]. In Eqs. (|28l) and (|29|) . we decomposed this transformation into two factors, 
one affecting the core states, and the other affecting the other states. 

The "self-similarity" of the wave function, evident in Eqs. f53|) - suggests that the effect of vortex encircling 
on the states \A X ) is analogous to its effect on the ground states |m) themselves. Clearly, the analog of the operators 

7i in the first generation of near-core states are the operators z\ X \ Thus, one expects 

B i:j = 2Z { pzf (56) 

and the associated operator identity 

2 7i7j = {2i ll X l ){2t 1] X j ){2z[ 1) zf ) ) (57) 

when applied to any ground state. Indeed, the proof of this identity follows from the fact that Yi — (l/V2)(Xi — iZ\ ■ ) 
annihilates any ground states so that, when acting within the subspace of ground states, one has = = 
1/2 - iXiZP or 2iX % zf ) = 1. 

An analogous picture emerges for the exchange of vortices. Here, we found that the explicit effect of exchange on 
the states |m, x) is equivalent to the action of the operator exp{nji'yj/2} exp{nXiXj/2}. On the other hand, the 
effect of the interchange on the ground states is given by the operator exp{7T7i7j/2}. Due to the "self-similarity" of 
the wave function, the connection between these two operators is expected to be furnished by 

B = exp{nZ^Z ( J 1) /2}. (58) 

Indeed, there exists a corresponding operator identity 

exp{7r 74 7 3 /2} = exp{7r 7i 7 i /2} exp{ 7 rX l X J /2} exp{^Z l (1) Z J (1) /2} (59) 

valid when acting within the subspace of ground states. Its proof uses the same ingredients as in the case of encircling. 

For both vortex encircling and exchange, this procedure can now be repeated for the near-core states obtained in 
higher iterations. In this way, it follows that the effect of vortex encircling or exchange on any core or near-core states 
is entirely contained in the geometric phases f2i. 

This discussion clarifies the difference between unitary transformations composed of zero-energy Majorana opera- 
tors, say, 72j-i72j) an d unitary transformations composed of non-zero energy Majorana operators, such as -^2j-i^2^ 

(k) (k) 

and ZfiLj.Ziy. When acting on a ground state, the first one leaves the system within the subspace of ground states, 
while the last two excite the system. Any unitary transformation that may be implemented by means of an adiabatic 
vortex braiding leaves the system in the ground state subspace. That is the case for 723-17231 whose effect on a 
ground state may be reduced to a series of phase changes associated with vortex braiding. By contrast, the effect of 
the other two operators may not be implemented by an adiabatic braiding of vortices. 

The number of generations L appearing in Eq. (|55|l is chosen such that the states (r) generated by our iterative 
process do not become extended enough to overlap. If the iterative process is carried out further, states from different 
vortices start overlapping, and therefore cannot be used as basis states in a single-particle basis. In principle, these 
states may be orthogonalized by means of a Gram-Schmidt-type process. The states resulting from that process, 
however, would not necessarily share the properties of the states (r) for k < L, namely their occupation may not 
be particle-hole symmetric, and the way they are affected by vortex braiding may be different. 

VII. SUMMARY 

The composite-boson theory of the fractional quantum Hall effect at filling fractions v = 1/m (m odd) employs 
a Chern-Simons transformation to map the electronic system to a superfluid of composite bosons. It explains the 
fractional statistics of the excitations of these states by describing them as charged vortices in that superfluid, and 
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the statistical phase as the effect of the charge of one vortex on the geometric phase accumulated by another vortex 
as it moves adiabatically. 

The Moore- Read theory describes the v = 5/2 state also as a superfluid. Again, charged excitations are vortices 
in this superfluid, and these vortices accumulate geometric phases as they move. However, in this case the effective 
bosons forming the superfluid are Cooper pairs of composite fermions. As a consequence, this state has excitation 
modes that involve the breaking of Cooper pairs into two composite fermions. These modes are solutions of the 
Bogolubov-deGennes equations. In the presence of well-separated vortices, each vortex gives rise to one zero-energy 
solution. These zero-energy solutions represent Majorana operators, operating on the core state of their vortex. For 
2A^ vortices, these zero modes lead to a 2 w -fold degeneracy of the ground state in the presence of vortices. The 2N 
core states define a 2 2Ar -dimensional Fock space, since each core state may be occupied or empty. We found that each 
ground state has an equal probability for each of the 2 2N possibilities to occupy the core states, and that the way 
these possibilities are superposed is such that (at least for ground states with a definite parity of the total particle 
number) the occupation of core states in one vortex is entangled with the occupations of all core states in all other 
vortices. We showed that it is the relative phases between the amplitudes for different occupations that distinguishes 
between the various ground states, and that it is these relative phases that vary as the vortices move. This allows the 
system to transform between ground states due to vortex braiding. It is interesting to observe that unlike the case 
of ordinary quantized Hall states, there seems to be no direct relation between the charge carried by the vortices and 
their non-abelian statistics. 

While the zero-energy Majorana modes affect only the occupations of the core states of the vortices, we showed that 
it is possible to extract some general information regarding the occupations of other single-particle states without 
explicitly solving the Bogolubov-deGennes equations. In particular, we constructed a set of single-particle states 
which are localized near the vortex cores and whose occupations are also maximally uncertain (all possibilities of 
their occupations are equally probable) and maximally entangled. The definition of these near-core states allowed us 
to explicitly relate our picture of nonabelian statistics, based on entanglement and geometric phases, to previously 
existing approaches based on representations of vortex braiding in terms of unitary transformations that create and 
annihilate particles in the core states. The elucidation of the role of quantum entanglement and geometric phases 
may be useful for understanding the effect of decoherence on non-abelian statistics, a question which is of relevance 
in the context of topological quantum computational 

The Moore-Read state is only the simplest example of proposed FQHE wave functions- 2 whose excitations satisfy 
nonabelian statistics. As in this example, nonabelian statistics is generally associated with a degeneracy of the ground 
state. However, the other proposed states lack a description in terms of second quantization, and in the absence of 
such a description it is hard to generalize our analysis for the Moore- Read state to the entire set of nonabelian states. 
We believe, however, that a second-quantization formulation of other nonabelian states would lead to a similar picture. 
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APPENDIX A: UNITARY TRANSFORMATIONS FOR VORTEX ENCIRCLING AND INTERCHANGE 

In this appendix, we review for completeness how to construct the unitary transformations U associated with 
particle exchange and encircling^ For this purpose, we analyze the time evolution of the many-particle state within 
the degenerate subspace of ground states as the vortices adiabatically traverse trajectories that start and end in the 
same set of positions. The unitary transformation is defined by the relation between the final state \ip(t — T)} and 
the initial state \ip(t — 0)), 
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m = T)) = um = o)). 



(Al) 



Correspondingly, the time evolution of the operators 7j spanning the degenerate subspace is given by 



li (t = T) = U li {t = Q)U^. 



(A2) 
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The time evolution of the operators ji can be readily read off from their definition Eq. ijSJ) and the explicit form Eq. 
(0) of the zero-energy spinors of the Bogolubov-deGennes equations. This can be used to identify the operators U up 
to an abelian phase. We note that with our choice of phase for the zero-energy spinors in Eq. JJJ, their Berry phase 
during adiabatic phase evolution vanishes and the phase evolution is given solely by the explicit monodromy. 

Before constructing U for exchange trajectories or encircling, we first construct unitary operators from the set of 
7i's which are helpful in giving explicit expression for the operators U. 



• The unitary transformation u k = V^lk adds a minus sign to all operators 7j with j ^ k. If j — k, it leaves the 
operator unchanged: 

t _ / Ik when j = k , . . 

klj k S otherwise ^ ' 

Obviously {u k ,u k >} = 25 kk >. 

• The unitary transformation uy — ji + jj interchanges 7i with "fj and adds a minus sign to all other operators 
7 fe with k ^ 

!ji when k = j 
7j when k = i (A4) 
— 7fc otherwise 



1. Winding trajectories 

We start with the case that all vortices move along closed trajectories. If m k is the number of windings of the fcth 
vortex around other vortices, then fl k changes by 1-Km k . In view of Eqs. and JSJ), this implies that the Majorana 
operator 7^ is multiplied by (— l) mfc . 

In the simplest case vortex 1 encircles vortex 2, leading to both 71 and 72 being multiplied by —1, with all other 
operators unchanged. This is a consequence of the fact that a phase factor 2-k of the order parameter leads to a phase 
7r (and thus a minus sign) for the fermionic operators. 

In the general case, we need to find a unitary transformation J7 w i n( j by requiring 

t _ J -7 fe when m k is odd . 
iv W ind7^wind ~\ lk when m k is even ' [A0> 

This is satisfied by the operator 

27V 

?7wiad=n U r> ( A6 ) 
fc=l 

where the operator ordering in C/ W i nc i is not important up to an overall minus sign. 



2. Exchange trajectories 

Exchange trajectories in which some of the vortices trade places are more complicated since the phase changes of 
fife associated with a particular trajectory do not only depend on the winding numbers, but also on the details of the 
trajectory and on the precise definition of the cut of the function arg(r) where its value jumps by 2tt. 

The simplest example is the interchange of two vortices. Inevitably one of the vortices crosses the cut line of the 
other vortex. Thus such an interchange involves an operator U12 that exchanges the positions of the two vortices and 
an operator u\ or ui that multiplies the appropriate vortex by —1. As a result, the unitary transformations generated 
by these exchanges are 

U 1 - 2 = _(l + U2Ul ) (A7) 
v2 

U 1 - 2 = _(l + Ul «2). (A8) 
v2 
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The transformation U 1 ^ 2 transforms ci — > C2 and c-i — > — ci while t/ 1 ^ 2 transforms c\ — ► —c^ and C2 — > ci. 

The combination of the motion of the first vortex from R4 to R2 and the motion of the second vortex from R 2 to 
Ri generates a closed curve. When that curve encloses a third vortex, the phase of that vortex changes by 27r, and so 
does the phase of one of the first two vortices. Altogether, this gives three phase shifts of 27r, and the transformations 
U 1 ^ 2 and U 1 ^ 2 have to be multiplied by either U3U1 or 113112- There are two possible outcomes to that. The first 
is U3 +U2), a transformation in which the phase of each of the three vortices is shifted by 2ir. The second is 
U3 (ui — 112), in which the phase of the third vortex is shifted by 2tt, and the other two 2-7T shifts are both given to 
one vortex. 



APPENDIX B: ORTHOGONALITY OF VORTEX STATES 



In this appendix, we discuss the scheme to generate the functions (r) in more detail and prove the orthogonality 
relation Eq. (J25J. We start by proving the orthogonality of the states u^°'(r) and Wj(r). In spinor notation, we can 
write 



in terms of the expansion coefficients 



(Bl) 



r 3 — 



1_ / U E 

V2\v E 



(0) 

„(°)l 



These coefficients satisfy the relation 



C 3 -e = [C 3 E Y 



(B2) 



(B3) 



whose proof uses the identity ue(i") = [i?_s(r)]*. The states io| (r) and w^ir) are obviously orthogonal due to their 
localization properties for i ^ j. For i = j, the orthogonality follows from 



Re(^ (1) l^ (0) ) 



w 



mi. 
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(B4) 



and 



T / (!)| (°)\ 

lm(^ \wl ) 







■ (0) \ 




V K (1) ]* 





(B5) 



Here, Eq. (|B4|I uses that by construction, the spinor zj 1 ' has zero overlap with the zero-energy spinors. Eq. I|B5|I uses 
the symmetry l|B3(l of the expansion coefficients C E . 

We now turn to the higher iterations (r) and proceed to prove the orthogonality Eq. (|42(l . Using the completeness 
of the cigcnfunctions of the BdG equations, we obtain the useful result 



w 



(0) 

j 

(0)i, 
3 ' 



(0) 
'V- 
3 



(B6) 



with S a hermitian operator, S — . 

For different vortices, i ^ j, the orthogonality follows again from the locality properties of the iu's. Thus, we focus 
on states from the same vortex, i = j, and drop the subscript labelling the vortex in the remainder of this appendix. 
In each iterative step k of the construction of the ui's, we need to prove that 



(w il) \w {k) ) =0 



(B7) 



for all / < k. Since this amounts to a proof by induction, we may exploit that two states w with indices smaller than 
k are orthogonal. 
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Starting with I — 0, we have 
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since by construction, the spinor has zero overlap with zero-energy spinors. Furthermore, 



For I = 1,2. 
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, k — 2, we obtain 

2Re(w (fc V (0 > 



w (l) 



w (k) 
w (k-i) 

1 / w 
~2\ M*" 1 ) 



(fe-i) 



«,(*-!) 
[«;(*-!)]• 



w (0) 



(-5o- a )t(-Str,) 



w (i-i) 
[wV-V]* 



(0) 



it) 



,('-!) 
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= — z 
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(B8) 



(B9) 



(BIO) 



(Bll) 



Finally, we note for I = k — 1 that Ke(w^\w^ k ^) = can be obtained in complete analogy to Eq. IjBlOfl and 
Im(t(/ fc ) |t// fe-1 ') = in analogy with Eq. (|B5(I . This completes the iterative construction of the functions w 
including the proof of the orthogonality relations Eq. H42[l . 
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